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Abstract 

For a decade and a half GPS Common-View time transfer has greatly served the needs of primary 
timing laboratories for regular intercomparisons of remote atomic clocks . However, GPS as a one- 
way technique has natural limits and may not meet all challenges of the comparison of the coming 
new generations of atomic clocks . Two-Way Satellite Time and Frequency Transfer (TWSTFT) is 
a promising technique which may successfully complement GPS, For hvo years , regular TWSTFT’s 
have been performed between eight laboratories situated in both Europe and North America , using 
INTELSAT satellites. This has enabled an extensive direct comparison to be made between these two 
high performance time-transfer methods. The performance of the TWSTFT and GPS Common-View 
methods are compared over a number of time-transfer links . These links use a variety of time-transfer 
hardware and atomic clocks and have baselines of substantially different lengths. The relative merits 
of the two time-transfer systems are discussed. 
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INTRODUCTION 


The performance of atomic clocks maintained at primary timing laboratories have improved 
considerably in recent years. There is now a challenge to develop suitable time and frequency 
transfer methods to exploit this improved performance. The standard method of intercomparing 
clocks contributing to International Atomic Time (TAI) is by common-view of Glob al Pos itioning 
System (GPS) satellites! 1 !. TWo-Way Satellite Time and Frequency Transfer (TWSTFT) has in 
recent years been developed as an alternative time and frequency transfer method! 2 !. TWSTFT 
as a two-way time-transfer method, offers many potential advantages over the existing one-way 
methods. Due to the symmetrical nature of the TWSTFT method, several sources of systematic 
errors are either eliminated or greatly reduced. These include errors associated with Earth 
station and satellite positions along with ionospheric and tropospheric delay errors. The use of 
directional antennas and high frequency transmissions enables low power transmissions to be 
made with relatively high carrier-to-noise ratios, resulting in high precision measurements. The 
downside is that TWSTFT instrumentation is somewhat more expensive. Satellite time must 
also be purchased on a commercial geostationary satellite. 

In this paper, a detailed study is presented of the comparison between regular 
TWSTFT and GPS common-view measurements. Measurements have been included from five 
European laboratories recorded over a period of two years. The (TWSTFT-GPS) differences 
obtained from each link were examined. Values of a t were calculated from the TWSTFT 
and GPS time transfers and also from the (TWSTFT-GPS) differences. Comparisons were 
made against a y values calculated from co-located atomic clock comparisons performed at NPL. 
Finally, discrepancies between the TWSTFT and GPS time transfers are explained in terms of 
changes in both instrumentation and environmental conditions at each laboratory. 


METHOD 

The international TWSTFT field trial experiment has been performed during the last two years 
using an INTELSAT satellite at 307°E! 3 1. Six European and two North American laboratories 
have been participating in this experiment. The instrumentation used at each location has 
been described previously I** 5 !. The results presented in this paper have been obtained from 
interoomparisons of data from five of the European Laboratories. These were the Technical 
University Graz, Graz, Austria (TUG), National Physical Laboratory, Teddington, UK (NPL), 
Van Swinden Laboratorium, Delft, the Netherlands (VSL), Forschungs- und Technologiezentrum, 
Deutsche Telekom, Darmstadt, Germany (FTZ), and Physikalisch-Technische Bundesanstalt, 
Braunschweig, Germany (PTB). The atomic clocks and TWSTFT Earth station instrumentation 
used at each location are summarized in Thble 1. The provision of cost-free time by INTELSAT 
enabled the TWSTFT measurements to take place. Initially the INTELSAT (VA-F13) satellite 
at 307°E was used, but this was replaced by the INTELSAT (VII-F6) satellite. A schedule 
of TWSTFT measurements has been performed three times per week, with each individual 
tim e transfer lasting for five minutes. GPS measurements were made according to the BIPM 
International GPS common-view schedules 22, 23, 24, and 25. All GPS time receivers involved 
in this study are single-channel, C/A code, NBS-type. 


DATA ANALYSIS 

The TWSTFT and GPS data sets are fundamentally different. The TWSTFT data consist of 
spot measurements of five minutes duration made every two or three days. In contrast, the GPS 
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readings are obtained from the mean of a series of up to 60 thirteen-minute measurements, 
spread throughout two days. The noise within a five-minute TWSTFT time time transfer is 
substantially lower than the noise within a block of 13-minute GPS data. Underlying both sets 
of measurements are the variations of the atomic clocks. These short- and medium-term clock 
variations are small in the case of the hydrogen maser, but significantly larger in the case of 
commercial cesium clocks. 

TWSTFT measurements were made on Mondays, Wednesdays, and Fridays. This measurement 
schedule resulted in a regular but unevenly spaced data set. A simple algorithm has been 
developed to calculate a good approximation to cr y 1*1 under these conditions. This algorithm 
has been successfully implemented. 

Time transfers were computed using the GPS common-view method for the period MJD 
(49354-49950). During this period the Block II satellites were permanently subjected to 
Selective Availability (SA), so strict common-views were required to remove the effect of the 
SA clock dither. All common views retained for this study fulfilled the following conditions: 15 
s common-view tolerance, 765 s minimum duration of the track, 20° minimum elevation angle 
for satellites. The 15 s tolerance for common-views was necessitated by a fault in the NBS 
type receivers which begin observations 15 s later than scheduled. There were between 25 to 
40 common-view tracks per day fulfilling these conditions. Values of the common-views were 
computed for the midpoints of the tracks. The coordinates of the GPS ground antenna were 
expressed in the ITRF88 reference frame with uncertainty ranging from 10 cm to 30 cmril. 
The coordinates for NPL, VSL, and FTZ were newly determined and the GPS data were 
corrected in post-processing. The distances between European time laboratories range from 
a few hundred kilometers to about one thousand kilometers. During the GPS common-view 
time transfer the errors due to broadcast satellite ephemerides, ionospheric, and tropospheric 
delays, are reduced to the level of 1 ns or lower !®»®1. Therefore, there is no need to use 
post-processed precise ephemerides and measurements of ionosphere and troposphere I 1 ®!. For 
each link, a Vondrak smoothing was performed on the values UTC(Labl)-UTC(Lab2), which 
acts as a low-pass filter with a cut-off period ranging from 0.5 day to 2 days depending on the 
pair of laboratories! 11 !. Those cut-off periods have been chosen as being approximately the limit 
between the short time intervals, where the measurement noise is dominant, and the longer 
intervals where the clock noise prevails. Finally, the smoothed values were inteipolated for the 
occurrence of the TWSTFT measurements. The Vondrak smoothing method is illustrated in 
Figure 1. The “cloud" of GPS data points are shown, along with the cuive resulting from the 
smoothing. The GPS links were differentially calibrated with an uncertainty of about 2 ns! 12 !. 


RESULTS 

Curves of the (PTB-NPL) time transfer made over a two-year period are shown in Figure 2. 
The offset between the two curves is due in part to the delay asymmetries of the TWSTFT 
instrumentation not being calibrated, and in part to an offset of 150 ns being added to clearly 
separate the two curves. Curves of the (TWSTFT-GPS) differences are shown in Figures 3, 4, 
and 5. Values of a y calculated for both the TWSTFT and GPS Common-View time transfers 
and (TWSTFT-GPS) differences are shown in Table 2. a y values were calculated with averaging 
times (r) of 2.3, 4.7, and 7 days. 

Several trends emerged. There is good agreement in the shape of the time-transfer curves 
obtained using the TWSTFT and GPS common-view methods. Values of the standard deviation 
calculated from the (TWSTFT-GPS) differences are shown in Table 3, both for the complete 
data set and for a sub-section. Outlying points that deviated substantially from the mean value 
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were removed before calculating the standard deviation using the points shown in Figures 3-5. 
The NPL-TUG differences exhibited the lowest standard deviation when calculated over the 
whole period. In order to compare (TWSTFT-GPS) differences over shorter periods, standard 
deviations have been calculated from subsets of the data which are free from rapid delay 
changes. The results obtained were encouraging. Standard deviations of between 1.4 ns and 
2.7 ns were obtained over periods ranging from 70 to over 300 days. The most stable operation 
occurred on the (PTB-NPL) link, where a standard deviation of 1.4 ns was obtained for the 
(TWSTFT-GPS) differences over a period of 350 days. 

Values of <r„ were calculated from TWSTFT, GPS, and (TWSTFT-GPS) data sets that contained 
only data collected on the same MJDs. When data were missing from one data set, the 
corresponding data were removed from the other data sets before processing. Any discrete 
delay steps occurring due to known instrumentation changes were removed before the <x y values 
were calculated. Large discrete delay steps of amplitude 20 ns and 30 ns were removed from 
the FTZ data sets before calculating a y . These steps occurred only occasionally within a data 
set and were not typical of the data scatter. Values of a y varied considerably between the 
TWSTFT links. The (PTB-NPL) time transfer was the most stable link. These results were 
attributed to the use of an active Sigma Tau hydrogen maser at NPL and the Primary Cesium 
clock (CS2) at PTB compared with the use of commercial HP5071A cesium clocks at the other 
laboratories. With averaging times (r) of 2.3 days or longer, the principal instability contributing 
to the <r y values was clock noise. This is explained below. 

For a given time transfer, values of a y were in almost all examples lower for the GPS common- 
view measurements when compared against the TWSTFT measurements. In most cases, the 
difference in a y values was quite small, but clearly significant. This difference was particularly 
noticeable on the most stable (PTB-NPL) link, with an averaging time (r) of 2.3 days. In almost 
all examples, the values of a y obtained from the (TWSTFT-GPS) difference were significantly 
lower than the a y values obtained from the individual time transfers. This again indicated 
that the major contribution to the time transfer a y values is from clock noise. A significant 
proportion of this noise cancels in the (TWSTFT-GPS) differences, due to the partial elimination 
of noise from the clocks. 

Despite the lower a y values obtained from the GPS time transfers, the conclusion should not 
be drawn that the GPS common-view method offers the best technique for clock comparison. 
The lower <r y values may be due to the choice of TWSTFT and GPS measurement schedules, 
rather than to an intrinsically higher accuracy of the GPS method. The a y values obtained 
from a TWSTFT are calculated from “spot" five-minute readings, made either two or three 
days apart. In contrast, the a y values obtained from a GPS common-view time transfer are 
calculated from “weighted means" of up to two days’ data, with approximately thirty satellite 
readings contributing to each day’s data. a y values calculated from these mean values will be 
lower even in the case where two perfect time-transfer systems are used. 

To illustrate the above point further, co-located clock comparisons have been made between two 
HP5071A commercial cesium clocks and an active hydrogen maser at NPL. One hundred days 
of measurements have been examined. Values of a y obtained from the comparisons are shown 
in Table 4, using averaging times ( t ) of 2, 4, and 8 days. The a y values were calculated from 
single readings, from the mean of 48-readings taken over the two days, and from the (single 
reading - mean readings) differences. The results show the advantage of taking measurements 
throughout the 48-hour period. With a two-day averaging time, the a y values calculated from the 
mean readings were substantially lower than the a y values calculated from the single readings. 
The values of a y obtained from the co-located measurements were comparable with, and only 
slightly lower than, the a y values obtained from the TWSTFT and GPS time transfers. These 
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results suggest that a large fraction of each time transfer a y value is due to clock noise. These 
results also show that the differences between the a y values obtained from the TWSTFT and 
GPS common-view systems is most likely to be due to the choice of measurement schedule, 
rather than any intrinsically better delay stability of the GPS system. 

Plots of (TWSTFT-GPS) differences show several trends. There are delay steps of several 
nanoseconds occurring within some of the instrumentation. Several, but not all, of these 
delay changes are associated with known delay changes of either the GPS or the TWSTFT 
instrumentation due to hardware replacement. Periodic delay changes with an annual period 
occur within the (VSL-NPL) differences. Temperature-dependent delay changes have been 
observed in previous studies of (TWSTFT-GPS) differences I 13 l. 


DISCUSSION 

There are several possible improvements that may be made to both the TWSTFT and GPS 
common-view time-transfer systems used in this experiment. Neither the 
TWSTFT system nor the GPS system are presently operating using optimum hardware. The 
GPS system may benefit from the use of dual-frequency multichannel receivers, which make 
optimum use of the available GPS signals. Work has already been performed to improve the use 
of the GPS system for time transfer by, for example, better satellite ephemeris determination, 
improved Earth station coordinate determination, ionospheric measurement, and tropospheric 
modelling. The main limitation to the performance of the common-view GPS method is the 
delay stability of the receiver instrumentation. The TWSTFT system may benefit from the 
use of more recently designed modems. Further improvements may be obtained from the 
optimization of the Earth station instrumentation to minimize the delay instabilities. Satellite 
simulators may be used to measure the Earth station delay asymmetries during a TWSTFT 
measurement session. 

The values of a y obtained using both the TWSTFT and GPS systems have been limited by the 
performance of the cesium atomic clocks at most locations. A parallel TWSTFT experiment 
has been taking place between Europe and North America. The combination of longer baseline 
time transfers and the possibility of operating with active hydrogen masers at both locations 
should make the study of these links of considerable interest. It will be of particular interest 
to examine the effects of ionospheric corrections, and precise ephemeris corrections applied to 
the GPS measurements, made over these relatively long links. 


CONCLUSIONS 

TWSTFT and GPS common-view methods have been shown to be capable of providing high- 
precision time transfers. Values of a y (r = 2.3 days) as low as 1.3 x io -14 and 1.8 x ltr 14 have been 
reported for GPS common-view and TWSTFT respectively. This difference in a y values has 
been attributed to the behavior of the clocks when interrogated using the different measurement 
periods used by the two systems. Periodic delay changes of period one year were also observed. 
These changes correlated with outdoor temperature variations. Further work is required to 
obtain the optimum performance from both systems. 
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Table 1 Atomic clocks and TWSTFT instrumentation. 


Laboratory 

TUG 

NPL 

VSL 

FTZ 

PTB 

Clock 

HP5071A 

Hydrogen 

Maser 

HP5071A 

HP5071A 

CS2 Primary 
Caesium 

TWSTFT 

Antenna 

1.8m 

2.4m 

3.0m 

1.8m 

1.8m 

TWSTFT 

Modem 

MITREX 

2500 

MITREX 

2500 

MITREX 

2500 

MITREX 

2500A 

MITREX 

2500A 


Table 2 Values of o y calculated from TWSTFT and GPS time transfer and (TWSTFT-GPS) 
differences. 


LINK 

TWSTFT Time 
Transfer o y xlO 14 

GPS Time Transfer 
o y xl0 14 

(TWSTFT - GPS) 
Difference o y xl0 14 


m 

t=4.7 

days 

t=7.0 

days 

mm 

t=4.7 

days 

t=7.0 

days 

t=2.3 

days 

t=4.7 

days 

t=7.0 

days 

NPL-TUG 

2.4 

b 

1.2 

2.3 

B 

1.2 

1.3 

0.6 

0.5 

VSL-NPL 


5.5 

4.5 

4.8 

5.3 

4.4 

2.8 


0.9 

VSL-TUG 

5.8 

4.5 

4.8 

5.6 

4.0 

4.4 

2.0 

B 

1.0 

FTZ-VSL 


5.3 

B 

4.4 

4.9 

4.6 

4.6 

B 

B 

FTZ-NPL 

3.6 

B 

2.5 

3.0 

2.5 

2.4 

2.0 



FTZ-TUG 

4.3 

2.7 

2.6 

4.0 

2.7 

2.5 

2.3 

B 

B 

PTB-FTZ 

3.8 

B 

2.4 

3.6 

2.4 

2.3 

0.9 

0.7 

B 

PTB- VSL 

4.3 

3.6 

3.4 

4.2 

wnm 

3.2 

2.6 

B 

1.0 

PTB-NPL 

1.8 

1.0 

1.0 

1.3 

0.8 

0.9 

1.7 

0.8 

0.6 

PTB-TUG 

2.9 

1.6 


2.6 

1.6 

1.4 

1.5 

0.9 

0.8 
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Table 3 Standard deviations calculated from the (TWSTFT-GPS Common-View) 

differences. 


LINK 

Standard 
Deviation 
(Complete Period) 

Duration of 
subset 
(MJDs) 

Standard 

Deviation 

(subset) 

NPL-TUG 

2.4 ns 

49471-49840 

1.5 ns 

VSL-NPL 

4.1 ns 

49707-49805 

1.4 ns 

VSL-TUG 

4.6 ns 

49590-49709 

2.7 ns 

FTZ-VSL 

14.5 ns 

49670-49754 

2.9 ns 

FTZ-NPL 

7.4 ns 

49670-49805 

1.6 ns 

FTZ-TUG 

9.1 ns 

49635-49805 

1.5 ns 

PTB-FTZ 

11.7 ns 

49657-49805 

1.6 ns 

PTB-VSL 

4.3 ns 

49670-49840 

2.4 ns 

PTB-NPL 

2.5 ns 

49600-49950 

1.4 ns 

PTB-TUG 

3.4 ns 

49567-49950 

2.3 ns 


Table 4 o y values calculated from co-located measurements made at NPL 



Single Readings 
o y xlO 14 

Mean Readings 

OyXlO H 

Difference 

OyXlO U 


t=2.0 

days 

t=4.0 

days 

t=8.0 

days 

t=2.0 

days 

t=4.0 

days 

t=8.0 

days 

x=2.0 

days 

r=4.0 

days 

m 

Maser-123 

1.8 

1.4 

0.7 

1.3 

1.3 

0.6 

1.7 

0.4 

0.3 

Maser-404 

3.2 

1.8 

0.9 

2.2 

1.4 

1.0 

1.6 

0.8 

0.3 

123-404 

3.6 

2.3 

1.0 

2.6 

1.9 

1.0 

1.8 

0.9 

0.3 


Maser = Sigma Tau Hydrogen Maser 
123 = HP5071A High performance clock 
404 = HP5071A Standard clock. 
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GPS DATA ANALYSIS 


Figure 1 
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lure 2 


TWSTFT AND GPS COMMON- VIEW TIME TRANSFERS 

PTB-NPL 


-2100 
-2200 
-2300 
-2400 
-2500 
-2600 
-2700 
-2800 

49300 49400 49500 49600 49700 49800 49900 50000 

MJD 



. _ GPS Time Transfer ^ TWSTFT + 150ns 


356 





Figure 3 


(TWSTFT - GPS COMMON- VIEW) DIFFERENCES 
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